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1 Introduction 



The canonical symmetry of an integrable system is a special discrete symmetry transforma- 
tion having a number of remarkable properties In particular, as we have shown in 
Ref. [Q, this symmetry, considered as a transformation of the phase space of the system, 
was a canonical transformation. In the same paper, considering the case of the nonlinear 
Schrodinger equation, we have established that the densities of the lowest conservation laws 
change under the action of the canonical symmetry by total space derivatives. Here we show 
that in this case all the conservation laws possess this property. From this it follows, in 
particular, that the whole hierarchy of the nonlinear Schrodinger equations is invariant with 
respect to the action of the canonical symmetry. 

Since, from our point of view, discrete transformations of the nonlinear differential equa- 
tions did not attract early the proper attention, we found it useful to give a short review of 
the basic facts on such transformations and their properties. This is done in sections 2 and 
3. Trying to make the paper self-contained but not too long, we decided to give here the 
proofs only for those facts that could not be found in the standard reference on the symmetry 
properties of the nonlinear differential equations |^. 

In Ref. to show that the canonical symmetry is a canonical transformation we used the 
method of generating functions. At the beginning of section 4 for the case of the nonlinear 
Schrodinger equation we give another proof of this fact based on the representation of the 
Poisson bracket with the help of the corresponding matrix differential operator The 
rest of section 4 is devoted to the proof of our main result concerning the behaviour of the 
densities of the conservation laws under the canonical symmetry. 

The summation over repeated indexes is implied. 



2 Hamiltonian Approach to Evolution Equations 

Here we recall the main facts on the Hamiltonian formulation of the nonlinear evolution 
equations We restrict ourselves to the case of two independent variables t, x and A 
dependent variables u"', a = 1, . . . , A. An evolution equation is a system of equations having 
the form 

u'} = K''{x,u,u,...,%), (2.1) 

where for any k m' denote a set formed by the fc-th derivatives u'" of the dependent variables 
over X. It is convenient to put = u"'. 

Consider the space A of the functions of the variables x, u"", and the derivatives of m° over 
X up to some finite order. For a general function of such type we use the following notation 

f[u]= f{x,u,u,...,%^). (2.2) 
The operator of the total derivative over x is defined as 

D^^ + t '-''4r- (2-3) 

fc=o du°- 
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Actually, when the operator D acts on an element of A we have only a finite number of 
members in series (p.3|) . For any two functions / and g we have 



D'^Ug) = E [^jD"^~\f)D\g). (2.4) 
The operator 

N(a) 

«N = E (^k[u]D\ (2.5) 

fc=0 

where G is called a differential operator of the order A^(a)0. Using Eq. ( |2.4| ), we can in 
an obvious way define the multiplication of differential operators. It can be shown that the 
multiplication of differential operators is an associative operation. 

The transposed (adjoint) operator of the differential operator a, given by Eq. (p.5|), is a 
differential operator a^, defined as 

^ E i-^fD'a, = E E (-1)™ D^-\a^)D\ (2.6) 

k=0 k=0 m=k \^/ 

Here and below we use the following convention: if the operator D, or some degree of the 
operator D, precedes some expression surrounded by brackets, it acts only on the expression 
in the brackets. 

A differential operator a is called symmetric (self-adjoint) if a"^ = a, and it is called 
skew-symmetric (skew-adjoint) if = —a. 

A matrix with the matrix elements being differential operators is called a matrix differ- 
ential operator. Any matrix differential operator A can be written as 

N{A) 

A=J2 ^kD\ (2.7) 

k=0 

where Ak are matrices with the matrix elements being elements of A. The transposed 
operator A^ of the matrix differential operator A, given by Eq. (|2.7| ), is defined as 



N{A) 

A^=Y: {-ifD'Al, (2. 

fc=0 

where is the transposed matrix of the matrix Ak. 

Introduce on A the operators of the variational derivatives over defined as 



fe=0 

It can be shown that 



= V(-1)'=Z^'^4-. (2.9) 
to c^- 



1^ = 0, a = l,...,A, (2.10) 



^The operators, involving negative degrees of I?, can also be defined, but we will not use such operators 
in this paper. 
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if and only if there exists a function g, such that 

f = Dg. (2.11) 

Hence, defining 

Ker 6/6u = {f eA\ Sf/du" = 0, a = 1, . . . , A} , (2.12) 

we can write 

Ker 5/5u = DA. (2.13) 

Two functions / and g from A are said to be equivalent if there exists a function h such 
that 

f = g + Dh. (2.14) 

It is easy to prove that we have defined an equivalence relation. The corresponding equiv- 
alence classes are called functionals. We denote the equivalence class containing a function 

/ by F. In accordance with the established tradition we write in this case F = J f dx. 

The space of all functionals is denoted by JF. If a functional F is given we denote by / an 
arbitrary representative from the corresponding equivalence class and call / the density of 
the functional F. 

A bilinear map {■,■}: x JF ^ JF, is said to be a Poisson bracket on JF if it satisfies the 
relations 

{F,G} = -{G,F}, (2.15) 
{F, {G, H}} + {G, {H, F}} + {H, {F, G}} = (2.16) 



for any functionals F, G and H. Equality ( p^.l6[ ) is called the Jacobi identity. 

A matrix differential operator J with the matrix elements J"'', a,b = 1, A, is said to 
be a Hamiltonian operator if the equality 

defines a Poisson bracket on JF. It is clear that the right-hand side of equality (|2.17|) does 
not depend of the choice of / and g. The Poisson bracket ( |2.16| ) can be shown to satisfy 
Eq. (|2.15| ) if and only if the operator J is skew-symmetric. The Jacobi identity imposes 
more complicated conditions on the operator J. We do not discuss them here, and only note 
that if J is a constant matrix then the Jacobi identity is valid. 

We say that evolution equation ( p.l|) can be written in the Hamiltonian form if there 
exists a Hamiltonian operator J and a function h E A such that 

j^a ^ (2.18) 

If this is the case, we write the equations of motion in the form 

ut = r'[u]^-^ (2.19) 
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and call equations (|2.19 ) the Hamilton equations. Actually Hamilton equations ( p. 19 ) are 



defined not by the function h but by the corresponding functional H, called the Hamiltonian 
of the system. 

A conservation law of evolution equation (|2.1|) is a relation of the form 

Pt = Df, (2.20) 



that should be valid in virtue of equation In Eq. (|2.2CI| ) p and / are elements of A 



which are called the density and the flow of the conservation law, respectively. If we have 
conservation law ( p.20| ) we say that the functional P is a conserved quantity, or an integral of 
motion. For the case of Hamilton equations (|2.19|) , the functional P, satisfying the equalityf] 



^ + {P,H} = 0, (2.21) 



is a conserved quantity. 



3 Differential Transformations 

Consider a transformation of the space described by x, dependent variables m", and their 
^-derivatives, that is determined by the relations 



'"^^ (3.1) 



We call such a transformation a differential transformation. If there exists a transformation 

u'' = ^^[u], (3.2) 

such that 



(p'^Mu]] = n^ (3.3) 

ve an invertible transforms 

defined by 



then we have an invertible transformation. The transformation law for the variables ^u"' is 



= =Z}>'^M. (3.4) 

Define the operator (p* that acts on a function f[u\ = f{x, u, u, . . . , 'u') according to the 
rule 

V*f[u] = /bN] = /(a;, ^[u], <p[u], ...,%[u]). (3.5) 
It is not difficult to get convinced that for any / G ^ we have 

Dy,*f = ip*Df. (3.6) 

Let X ^ -^^j i- 6. X is a set of functions ^ ^ = 1, • • • , The Frechet derivative of x 
is the I X y4-matrix differential operator x', defined as H] 



^\[uY[u]^^X> + ev[u]] 



(3.7) 

e=0 



^We allow P to depend explicitly on t. 
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for any v"- ^ A, a = 1, . . . , A. From this definition we get 



x\ = E §rD'- (3.8) 



In particular, for f&Awe have 



It follows from this equality that 



oo or 

= E (3-9) 



1^ = (3-10) 
Let X £ -^^j denote the set of functions (p*x\ ^ = 1; • • • ? by ip*x- Let us show that 

{^*xy[u]=x'[M¥[u]. (3.11) 

Using Eq. ( ^TSf ) we get 

(V'-Xf. = f V (II;) gfl". (3.12) 

n,m=0 \0 U"/ OU"- 

It is not difficult to show that 

A-D- = th)D-^J:^ n>ra, (3.13) 

^ D- = t(i\D-'J-, n<m. (3.14) 



ou"- k=0 \'^/ o u 

From this equalities and Eq. ( p.4|) it follows that 



=0 \OU°/ \n=Ofc=0 n=m+lk=0 \ / \0 U / 



Now using the identity 



we can reduce Eq. ( p. 151 ) to 



m n mm 



EE = EE' (3-16) 

n=Ofc=0 k=On=k 



n,m=0 \d / k=Q / KOU"'/ 

^t^-i^y-t^D". (3.17) 

That was to be proved. 

If differential transformation (|3.1|) is invertible, then using Eqs. (|3.3|) and ( p.ll|) we get 

i^'\[^[u]y\[u] = 61 (3.18) 
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From Eq. ( 3.11|) for any / G ^ we also have 
Taking into account Eq. ( p.lO| ), we get 



(3.19) 



(3.20) 



Consider now the behaviour of Hamiltonian equations (|2.19|) under differential transfor- 
mations. Note first that from Eq. ( |2.19| ) it follows that 



6u^ ■ 



For the transformed solution u"" we have 



^<^"N{")6 fa r 1 Tfccr .i^^N 



~a "Sr^ '^r' l"'iy">b la r n jbcr i 

= —J^'^t = ^ b[u]J [u\ 

n=0 OU'' 



6u^ 



(3.21) 



(3.22) 



Let us suppose that we deal with an invertible transformation, then we can define the 
differential operator J and the function h by 



/i[(y9['u]] = h[u]. 



In this case from Eq. ( p.20|) we get 



6h[u] 6{p*h[u] ,j, 6h[ip[u]] 



Hence, equality (|3.22|) can be written as 



(3.23) 
(3.24) 



(3.25) 



.a Tab r-r/^M 



(3.26) 



Thus, if we have a solution of equations ( |2.19| ), then after transformation ( p.l| ) we get a 
solution of equations ( 3.26| ). 

In the case when J = J, i. e. when 



(3.27) 



we call transformation ( ^.1| ) a canonical transformation. If, in addition, h — h & Ker 6/6u, 
i.e. 

h[ip[u]] - h[u] e Ker 5/5u, (3.28) 

then we call transformation ( p.l|) a differential symmetry transformation. In this case equa- 
tions ( |3.26| ) have the same form as equations ( ^.19| ). 
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4 Nonlinear Schrodinger Equation 

In fact, here and Ref. we consider the following complex extension of the nonlinear 
Schrodinger equation: 

ig + g" - 2erg^ = 0, ir - r" + 2er'^q = 0, (4.1) 

where q and r are arbitrary complex functions of the variables x and t, e is the coupling 
constant. In Eq. (O) and below dot and prime mean the partial derivative over t and x, 
respectively. The canonical symmetry for this system has the form ^, 

q = — , r = er^q — r" H . (4-2) 

er r 

The inverse transformation is 

q = e(fr — q" + — , r = —. (4.3) 
q eq 

Denoting = q, u'^ = r, we can write equations ( |4.1| ) in Hamiltonian form ( |2.19| ), where 

.^(r;). (4.4, 

h[q, r] = r'q' + er'^q^. (4.5) 
The differential operator 99', corresponding to transformation ( |4.2|) , can be written as 

where 

« = 0, /3 = -A' 7 = er^ (4.7) 

12 
2 
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(5 = -D^ + 2-D + 2erg - ^. (4.^ 



It is easy to show that in this case equality ( |3.27|) is valid, hence transformation ( |4.2|) is a 
canonical transformation. 

In Ref. m we have shown that for /i, given by Eq. ( [4.6| ) we have 



% r ^ 



/i[g,r] - /i[g,r] = D ( irg - — — ) . (4.9) 



Thus, transformation (^4.2| ) is a symmetry transformation. 

The densities of the conservation laws for the system, described by Eq. ( [4.1| ), are given 
by the formula [0, |] 

Pn[g,r] = rw„[g,r], n=l,2, (4.10) 
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where 

Wi[q,r]=q, (4.11) 
and the quantities m„ for n > 1 can be found from the recursive relation 

n-1 

Wn+1 = -iD{Wn) + er ^ WmWn-m- (4.12) 

m=l 

Multiplying Eq. ( [4.12| ) by r and using Eq. ( [4.10|) , we get the following recursive relation for 
the densities 

ir' 

Pn+l = -iD{pn) H Pn + ^ X! PmPn-m, (4.13) 

m=l 

with the initial condition 

Pi[q, r] = rq. (4-14) 
It is convenient for our purposes to introduce a fictitious conservation law 

zr' 

Po[g,r] = — (4.15) 
er 

and write Eq. (|4.13| ) as 

n-1 

Pn+l = -iD{pn) + e XI PmPn-m- (4.16) 
m=0 

Consider now the behaviour of the densities p„ under transformation ( |4.2| ). Denote by 
/„ the difference between the transformed density and the initial one, i. e. 

/n[g,r] =p„[g, r] -p„[g,r]. (4.17) 

For n = 0, 1 we have 

e \ tD{po) +pi J e 
From Eq. ( [4.16| ) we get for n > 1 the following recursive relation 

n-1 

fn-\-l ^-D(/^n) ~l~ ^ ^ i^fmPn—m Pmfn—m fmfn—m) • (4.19) 
m=0 

Using this relation, we find subsequently 

/2 = ^D [p^ + ^pi^ , (4.20) 

/s = (^p2 + epopi + jpI^ , (4.21) 

/4 = tD (p3 + epm + ^p\ + e'PoPi + jPo ) • (4-22) 
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Thus, it is natural to suppose that for n > 1 

/„ = D{dn), 

where the quantities dn are given by 



n I 



m,=l 



E 



Pil • • • Pi„ 



ilH \-i^+m=n 



Let us prove this statement. 

From Eqs. (|]^ and it follows that 



(4.23) 



(4.24) 



/n = ^Eh™"' E I^(PnK---R 

m=l \ iiH (-j„+m=n 



It is not difficult to show that 



n-l 



fn+l = iD{pn) + e E Pmfn-m- 
m=0 



Hence, it is enough to show that the equality 



n-l 



iD{pn + fn) =eJ2 fniiP 



m=0 



is true. From Eq. ( [4.25| ), using repeatedly Eq. ( [4.16| ), we get 

Pn + fn= {iD{Po) +Pl)c„-1, 

where Cq = 1, and for n > 1 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



E 



E 



Ph---Pir. 



m=l \ iiH yi^J^m=n 



(4.29) 



Taking into account Eq. ( |4.29| ) and the explicit expression for /o, we reduce Eq. ( [4.27| ) to 
the relation 

n-l 

D{Cn-l) = -it ^ fmCn-l-ra- (4.30) 
■m=l 

The validity of Eq. (^1301) follows directly from Eqs. (|]2|) and iWM) - 

Thus, we have shown that in the case of the nonlinear Schrodinger equation the densities 
of the local conservation laws change under the action of the canonical symmetry by total 
space derivatives. We call such densities quasi-invariants of the differential transformation 
under consideration. A direct consequence of the fact, we have proved, is the invariance of 
the higher nonlinear Schrodinger equations under the action of the transformations, given 
by Eq. 
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5 Conclusion 



It is likely that only the densities of the local conservation laws of the nonlinear Schrodinger 
equation are quasi-invariant under the action of transformations ( |4.2| ). If it is the case, then 
it is interesting to find the conditions under which the invariants of some discrete differential 
transformation, being a canonical transformation, are involutive with respect to the Poisson 
bracket. 

We see, that the condition of the quasi-invariance under the canonical symmetry can be 
considered as a characteristic of the densities of the local conservation laws. Note that this 
fact does not give us a constructive method to find such densities. We can use for this the 
recursion operator, which can be constructed with the help of a relevant Hamiltonian pair 
IP, corresponding to the considered system of nonlinear equations. In the next paper we 
show, that the notion of the canonical symmetry is useful in looking for such operators. 
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